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ABSTRACT 



Aims. The goal of this work is to investigate the Faraday rotation measure (RM) of radio galaxies in hot galaxy clusters in order 
to establish a possible connection between the magnetic field strength and the gas temperature of the intracluster medium. 
Methods. We performed Very Large Array observations at 3.6 cm and 6 cm of two radio galaxies located in A401 and Ophiuchus, 
^ [ a radio galaxy in A2142, and a radio galaxy located in the background of A2065. All these galaxy clusters are characterized by 
high temperatures. 

Results. We obtained detailed RM images at an angular resolution of 3^^ for most of the observed radio galaxies. The RM 
images are patchy and reveal fine substructures of a few kpc in size. Under the assumption that the radio galaxies themselves 
have no effect on the measured RMs, these structures indicate that the intracluster magnetic fields fluctuate down to such 
^— ^ ■ small scales. These new data are compared with RM information present in the literature for cooler galaxy clusters. For 
^sj , a fixed projected distance from the cluster center, clusters with higher temperature show a higher dispersion of the RM 
i/^ ■ distributions {(trm)^ mostly because of the higher gas density in these clusters. Although the previously known relation 
between the clusters X-ray surface brightness [Sx) at the radio galaxy location and grm is confirmed, a possible connection 
between the (trm — Sx relation and the cluster temperature, if present, is very weak. Therefore, in view of the current data. 



o 



^-^ it is impossible to establish a strict link between the magnetic field strength and the gas temperature of the intracluster medium. 
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1. Introduction the so-called Faraday rotation measure {RM). Therefore, 

information on the intracluster magnetic fields can be ob- 

The knowledge of the magnetic field properties in cms- ^ • -i • • - • --lu v u ^- r ^i u -^ 

^ ^ ^ ^ tamed, m conjunction with X-ray observations ol the not 

ters of galaxies has significantly improved in recent years ,i i , i ^ • c ^\. -Di\/r c a- i • • 

^ ^ r / S^s, through the analysis ol the RM ol radio galaxies m 

thanks to different methods of analysis (see e.g. reviews by ^u u i ^ • -n. i ^ ^ ^\^ ^ 

_ ^ n ^ the background or m the galaxy clusters themselves. 
Govoni & Feretti 2004, Feretti & Giovannini 2007, Ferrari 

et al. 2008, and references therein) . The presence of a mag- The RM studies of radio galaxies have been carried out 

netized plasma between an observer and a radio source on both statistical samples (e.g. Lawler & Dennison 1982, 

changes the properties of the incoming polarized emission. Vallee et al. 1986, Kim et al. 1990, Kim et al. 1991, Clarke 

In particular, the position angle of the linearly polarized et al. 2001) and individual galaxy clusters and galaxy 

radiation rotates by an amount that is proportional to the groups by analyzing detailed high resolution RM images 

path integral of the magnetic fields along the line-of-sight (e.g. Perley & Taylor 1991, Taylor & Perley 1993, Feretti 

times the electron density of the intervening medium, i.e. et al. 1995, Feretti et al. 1999a, Feretti et al. 1999b, Govoni 

et al. 2001, Taylor et al. 2001, Eilek & Owen 2002, Govoni 
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et al. 2008, Guidetti et al. 2010, Bonafede et al. 2010). 
These data are usually consistent with central magnetic 
field strengths of a few /iG, but stronger fields, with val- 
ues exceeding :^10/iG, are derived in the inner regions of 
relaxed cooling core clusters. 

For both merging and relaxed clusters, the RM dis- 
tribution seen over extended radio galaxies is generally 
patchy, indicating that the intracluster magnetic fields 
are not regularly ordered, but instead they have turbu- 
lent structures on linear scales as low as 10 kpc or less. 

In the last years, important progress has been made 
in the analysis of the magnetic field structure in galaxy 
clusters. In particular, it has been established that the 
magnetic field is likely to fluctuate over a wide range of 
spatial scales and, in a few galaxy clusters containing radio 
sources with very detailed RM images, the magnetic field 
power spectrum has been estimated (e.g. EnBlin & Vogt 
2003, Vogt & EnBlin 2003, Murgia et al. 2004, Vogt & 
EnBlin 2005, Govoni et al. 2006, Guidetti et al. 2008, Laing 
et al. 2008, Guidetti et al. 2010, Bonafede et al. 2010). 

An increasing attention is given in the literature to 
the connection between the magnetic field strength and 
the cluster gas density and temperature (e.g. Kunz et al. 
2010). SPH simulations (Dolag et al. 1999, 2002, 2005a), 
indicated that the magnetic field intensity declines from 
the center of the cluster outward with a rough dependence 
on the thermal gas density. Such magnetic field profiles 
were confirmed by other simulations based on different 
numerical techniques (e.g. Briiggen et al. 2005, Dubois & 
Teyssier 2008, Colhns et al. 2009). SPH simulations also 
predicted a very steep correlation of the mean magnetic 
field with the cluster temperature, e.g. B ocT'^ (Dolag et 
al. 1999). However, there is also the possibility that the 
magnetic field strength is only mildly dependent on the 
cluster temperature as recently pointed out by Donnert 
et al. (2009). Kunz et al. (2010) predicted a dependence 
of the magnetic field strength on the cluster temperature 
in cool-core clusters {B oc Ue T^/^ in cool-core clusters, 
B (X Tie in isothermal clusters). 

Dolag et al. (2001) showed that the correlation between 
two observable parameters, the dispersion of the RM dis- 
tribution {(Jrm) and the cluster X-ray surface brightness 
in the source location (5'x), is expected to reflect a correla- 
tion between the cluster magnetic field and the gas density. 
Indeed, measuring the slope of the cf^m — Sx correlation 
it is possible to infer the trend of the magnetic field versus 
the gas density. In addition, as predicted by cosmological 
magneto- hydrodynamic simulations (Dolag et al. 1999), 
galaxy clusters should have different central magnetic field 
strengths depending on their dynamical state and hence 
temperature, leading to an offset of this correlation for hot 
clusters. 

In this work we present sensitive polarimetric Very 
Large Array (VLA) radio observations at multiple fre- 
quencies of polarized radio sources embedded in a set 
of hot, nearby clusters of galaxies. These data are com- 
pared with RM data taken from literature, in order to 



investigate a possible connection between the magnetic 
field strength and the gas temperature of the intraclus- 
ter medium. Moreover, a detailed analysis of the Faraday 
rotation of the most interesting sources presented in this 
work will allow us, in a forthcoming paper, to determine 
the cluster magnetic field properties, by following the ap- 
proach proposed in Murgia et al. (2004) and applied in 
the cases of A2255 (Govoni et al. 2006), A2382 (Guidetti 
et al. 2008), and Coma (Bonafede et al. 2010). 

The paper is organized as follows: in Sect. 2 we de- 
scribe the selection criteria of the targets, we illustrate 
the radio observations, and the data reduction. In Sect. 3 
we show the results of the total intensity and polarization 
properties of the radio galaxies. In Sect. 4 we show the 
RM images. In Sect. 5 we investigate the magnetic field - 
gas temperature connection by comparing these new RM 
data and some other high quality data taken from litera- 
ture with the cluster temperature. Finally, in Sect. 6 we 
draw the conclusions. 

Throughout this paper we assume a ACDM cosmology 
with Ho = 71 kms-^ Mpc-\ n^ = 0.27, and Qa = 0.73. 

2. VLA observations and data reduction 

We selected a small sample of radio galaxies embed- 
ded in a set of nearby, high temperature galaxy clusters. 
Starting from the X-ray fiux limited sample by Edge et al. 
(1990), we considered all the clusters with a temperature 
T > 8 keV, a redshift z < 0.1, and declination S > -40° 
(to be observed by the VLA). Among the resulting nine 
clusters (Ophiuchus, Coma, A2319, A754, A2142, A401, 
A644, A2065, A2255) we skipped Coma and A2255 be- 
cause they were already analyzed in other RM projects 
(Feretti et al. 1995, Govoni et al. 2006, Bonafede et al. 
2010), and A754 because of its complex X-ray morphology. 
Of the remaining clusters, we searched in the NRAO VLA 
Sky Survey (NVSS; Condon et al. 1998) for bright and ex- 
tended radio galaxies. We selected all the radio galaxies 
located within 15' from the clusters center showing in the 
NVSS an angular extension LAS > V and a fiux density 
at 1.4 GHz > 60 mJy. There were no radio galaxies in 
A644 that meet the above selection criteria. In A2319 a ra- 
dio galaxy was selected, but the data were corrupted and 
not useful for our analysis. Thus, the final list of the new 
observed clusters, together with their redshift, the angular 
to linear conversion, the ASCA cluster temperature, and 
the corresponding temperature reference is given in Table 
[H We note that although A2065 is reported by Edge et 
al. (1990) and David et al. (1993) with a high temper- 
ature, other X-ray investigation found for this cluster a 
temperature T < 8 keV (e.g. Markevitch 1998). 

Here we present the results of multifrequency, polari- 
metric observations of these sources obtained with the 
VLA. The details of the observations are provided in Table 
[21 The sources were observed at 4885 MHz and 4535 MHz 
(4835 MHz for Ophiuchus) within the 6 cm band and at 
8085 MHz and 8465 MHz within the 3.6 cm band. A401 
was observed within each band in both the B and C config- 
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Table 1. List of hot, nearby clusters of galaxies hosting the selected targets. 



Cluster 


z 


kpc/^^ 


T 

[keV] 


Ref. 


A401 


0.074 


1.39 


8.3 ±0.5 


Markevitch (1998) 


A2142 


0.091 


1.67 


8.8 ±0.6 


Markevitch (1998) 


A2065 


0.073 


1.37 


5.4 ±0.3 


Markevitch (1998) 


Ophiuchus 


0.028 


0.55 


10.26 ±0.32 


Fukazawa et al. (1998) 



Col. 1: Cluster Name; Col. 2: Redshift; Col. 3: Angular to linear conversion; 

Col. 4: Cluster temperature measured with ASCA; Col. 5: Temperature reference. 



Table 2. Details of the VLA observations. 



Source 


R.A. 


Decl. 


Frequency 


Bandwidth 


Conf. 


Time 


Date 




J2000 


J2000 


MHz 


MHz 




hours 




A401A 


02 58 32 


±13 34 16 


4885/4535 


50 


B 


1.2, 1.5 


2002-Jul-13, 2003-Dec-22 








4885/4535 


50 


C 


1.2 


2002-NOV-04 








8085/8465 


50 


B 


1.2 


2002- Jul- 13 








8085/8465 


50 


C 


1.2, 1.4 


2002-NOV-04, 2004- Apr- 16 


A401B 


02 59 15 


±13 27 36 


4885/4535 


50 


B 


1.2, 1.5 


2002-Jul-13, 2003-Dec-22 








4885/4535 


50 


C 


1.2 


2002-NOV-04 








8085/8465 


50 


B 


1.2 


2002-Jul-13 








8085/8465 


50 


C 


1.2, 1.4 


2002-NOV-04, 2004- Apr- 16 


A2142A 


15 58 13 


±27 16 24 


4885/4535 


50 


B 


3.3 


2003-Dec-23 








8085/8465 


50 


C 


3.1 


2004- Apr- 13 


A2065A 


15 22 46 


±27 55 26 


4885/4535 


50 


B 


3.3 


2003-Dec-23 








8085/8465 


50 


C 


3.1 


2004- Apr- 13 


OPHIB 


17 11 54 


-23 09 23 


4885/4835 


50 


BnA 


1.6 


2003-Oct-07 








8085/8465 


50 


CnB 


1.6 


2004-Feb-08 


OPHIC 


17 12 09 


-23 28 43 


4885/4835 


50 


BnA 


1.6 


2003-Oct-07 








8085/8465 


50 


CnB 


1.6 


2004-Feb-08 



Col. 1: Source; Col. 2,3: Pointing position (R.A., Decl.); Col. 4: Observing frequency; 
Col. 5: Observing Bandwidth; Col. 6: VLA configuration; 
Col. 7: Time on source; Col. 8: Dates of observation. 



urations. A2142 and A2065 were observed in the B config- 
uration within the 6 cm band and in the C configuration 
within the 3.6c?72 band. Given its southern declination, 
Ophiuchus was observed with the VLA in the hybrid con- 
figurations BnA and CnB at 6 cm and 3.6 ctti, respectively. 
All observations were made with a bandwidth of 50 MHz. 



Calibration and imaging were performed with the 
NRAO Astronomical Image Processing System (AIPS) 
following the standard procedure: Fourier- Transform, 
Clean and Restore. Self-calibration was applied to remove 
residual phase variations. In the case of A401 the {u^v) 
data at the same frequencies but from different configura- 
tions were first handled separately and then combined. 



Images of the Stokes parameters /, [/, and Q have 
been obtained for each frequency separately. Images of 
polarized intensity P = (Q^ + /7^)^/^ (corrected for the 
positive bias), fractional polarization FPOL = P/I and 
position angle of polarization ^ = 0.5 tan~^(/7/(5) were 
derived from the /, Q, and U images. 



3. Total intensity and polarization properties 

In the following we give a brief description of total inten- 
sity and polarization properties of the individual sources. 

For each galaxy cluster the radio contours obtained 
from the NVSS overlaid on the ROSAT HRI X-ray images 
are shown (Figs. 1, 4, 6, 8). The NVSS images are at 
1.4 GHz and have an angular resolution of 45'^ The X- 
ray images are in the 0.1 — 2.4 keF band and have been 
convolved with a Gaussian with a = 16'^ High resolution 
images of the radio galaxies, in the 6 cm band, overlaid 
on optical red plate are inset. Every radio galaxy shows 
an optical counterpart in the DSS4^, except for A2065A 
which is likely a cluster background source. 

In addition, for each source total intensity and polar- 
ization images at 3.6 cm and 6 cm, convolved with the 
same beam of 3'^ are shown (Figs. 2, 3, 5, 7, 9, 10). 
Contours represent total intensity while vectors represent 
the orientation of the projected E-field and are propor- 
tional in length to the fractional polarization. In the frac- 
tional polarization images we considered as valid pixels 
those where the fractional polarization FPOL was above 
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Sa^poL' Sample plots of the trend of the fractional polar- 
ization as a function of the observing frequencies, calcu- 
lated in different regions of the sources are also shown. 

We note that in some sources we found some unpolar- 
ized areas. This is likely due to a combination of sensitivity 
limits, small-scale tangled fields, and depolarization. The 
calculation of the average fractional polarization may be 
biased if these depolarized regions are ignored. To reduce 
any possible bias, we calculated the fractional polarization 
over the sources in all the pixels where the total intensity 
signal is above Saj at all frequencies. 

The relevant parameters of the /, Q, /7 images, con- 
volved with a beam of 3'^ are listed in Table [3l The total 
flux densities, are also reported. The total flux densities 
were calculated, after the primary beam correction, by in- 
tegrating the total intensity surface brightness, down to 
the noise level. 



3.1. Abell 401 

The ROSAT HRI X-ray image presented in Fig. [T] shows 
the X-ray emission in the 0.1 — 2.4 keV band of the cluster 
of galaxies A401. In the southwest, at a projected distance 
of about 3 Mpc, not shown in the field of view of the image, 
lies the A399 system. The X-ray excess and the slight tem- 
perature increase in the region between the two clusters in- 
dicate a physical link between this merging pair of clusters 
(e.g. Fujita et al. 1996, Fabian et al. 1997, Markevitch et 
al. 1998, Sakelliou & Ponmann 2004, Bourdin & Mazzotta 
2008). Recent XMM X-ray analyses (Sakelliou & Ponman 
2004, Bourdin & Mazzotta 2008) pointed out that neither 
of the two clusters contain a cooling core and their central 
regions are nearly isothermal, with some small-scale in- 
homogeneities. However, the reasonably relaxed morphol- 
ogy of the clusters and the absence of major temperature 
anomalies argue against models in which A399 and A401 
have already experienced a close encounter (see also Fujita 
et al. 2008a). 

A401 is known to contain a faint radio halo around 
the cD galaxy at the cluster center (Harris et al. 1980, 
Roland et al. 1981, Giovannini et al. 1999, Bacchi et al. 
2003). We recently found (Murgia et al. 2010a) that the 
central region of A399 is also permeated by a diffuse low- 
surface brightness radio emission which we classified as 
radio halo. Indeed, given their proximity to each other, 
the interacting pair of clusters of galaxies A401 and A399 
can be considered as the first example of double radio halo 
system. 

The periphery of A401 hosts a number of tailed ra- 
dio galaxies. The head-tail morphology is probably caused 
by a sweeping of the radio emitting material due to the 
relative motion between the gas and the radio galaxy. In 
this work we have analyzed those labeled with A401A and 
A401B in Fig. [H 

The distribution of spectroscopic redshifts near the 
pair of clusters A399 and A401 shows only one peak. This 
is because the velocity dispersion for individual clusters is 



larger than the apparent velocity separation between A401 
and A399 (Hih & Oegerle 1993, Oegerle & Hih 1994, Yuan 
et al. 2005). On the basis of their radial velocity and their 
projected distance from A401, the two galaxies associated 
with the radio sources A401A and A401B are classified as 
cluster members. 

A401A is located to the west of the cluster cen- 
ter at a projected distance of about 5.8^ The high 
resolution image at 6 cm reveals a peculiar morphol- 
ogy. A compact slightly distorted double structure of 
about 20'^ in size whose optical counterpart is located at 
R.A.(J2000)=02h58m32.0s, Decl.(J2000)=13°34'20'^ and 
a low surface brightness tail of more than 2' in size with 
a morphology which appears not directly connected to 
the compact emission. The total flux density of the radio 
galaxy, both at 6 cm and 3.6 cm, is indicated in Table 3. 
The compact component has a flux density of 134 ±4 mJy 
and 86 ± 3 mJy at 6 cm and 3.6 cm respectively leading 
to a mean spectral inde:>o (^ — 0.7±0.1, typical for radio 
galaxies. In the tail the spectral index is very steep {a 
:^ 1.6±0.1). In fact the flux density is 25 ± 1 mJy and 
9.0 ± 0.3 mJy at 6 cm and 3.6 cm respectively. 

The head tail A401B is located at a projected dis- 
tance of about 8.5' to the southeast of the cluster cen- 
ter with the head pointed versus the south of the cluster 
and the tail elongated versus the north. The total flux 
density of the radio galaxy, both at 6 cm and 3.6 cm, is 
indicated in Table 3. The high resolution image at 6 cm 
shows the core, in the position R.A.(J2000)=02h59ml4.9s, 
Decl.(J2000)=13°2ri2'', coincident with a galaxy. At 
6 cm. the source is about 1' in size. 

Figs. [2] and [3] show the polarization images at 8465 
MHz and 4535 MHz, with an angular resolution of 3^' x 3" . 

The polarization is clearly detected at both frequencies 
in the bright compact double structure of A401A. The tail 
of A401 A has a steep spectrum, consequently, the total in- 
tensity is drastically reduced at 8465 MHz and the polar- 
ized signal is detected only in the brightest regions of the 
tail. By considering all the pixels where the total intensity 
signal is above 5<ji at all frequencies, the compact dou- 
ble structure of A401A has a mean fractional polarization 
of 12 ± 3% both at 8465 and 4535 MHz. In the tail the 
fractional polarization is higher than in the compact com- 
ponent reaching a mean fractional polarization of 36 ± 5% 
both at 8465 and 4535 MHz. In the bottom panels of Fig. 
[2] we show the trend of the fractional polarization as a 
function of the observing frequencies, at different source 
locations. The fractional polarization does not decrease 
significantly in the analyzed range of wavelengths. 

A401B is polarized at both frequencies up to a dis- 
tance of about 30'^ The mean fractional polarization is 
20 ± 3% at 8465 and 18 ± 3% at 4535 MHz. Therefore, 
the depolarization is negligible. This is also confirmed in 
the bottom panels of Fig. [3] where we show the trend of 
the fraction polarization as a function of the observing 



S(z^)oc ly ^ with a=spectral index 
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Table 3. Information on total intensity and polarization images at 6 cm and 3.6 c?n 



Source 


A 

(cm) 


Beam 

n 


a(I)* 
(mJy/beam) 


(mJy/beam) 


a(U) 
(mJy/beam) 


Flux density 

(mJy) 


A401A 


3.6 


3.0x3.0 


0.021 


0.018 


0.018 


95.0±3.0 




6 






0.031 


0.027 


0.028 


159.0=b5.0 


A401B 


3.6 






0.025 


0.025 


0.024 


36.0±1.0 




6 






0.033 


0.028 


0.029 


55.0±2.0 


A2142A 


3.6 






0.017 


0.015 


0.016 


22.0±0.7 




6 






0.029 


0.030 


0.028 


37.0±1.0 


A2065A 


3.6 






0.018 


0.015 


0.015 


19.0db0.6 




6 






0.038 


0.030 


0.029 


32.0±1.0 


OPHIB 


3.6 






0.044 


0.038 


0.036 


46.0±2.0 


OPHIC 


6 
3.6 






0.115 
0.036 


0.095 
0.036 


0.086 
0.039 


58.0=b2.0 
13.0±0.4 




6 






0.089 


0.092 


0.093 


23.0±0.7 



Col. 1: Source; Col. 2: Observation wavelength; Col. 3: Beam; 
Col. 4, 5, 6: RMS noise of the /, Q, U images; Col. 7: Flux density. 
* Note that in this table we indicate the values of cr(I), cr(Q), cr(U) 
8465 MHz and 4535 MHz (4835 MHz for OPHIB and OPHIC) at 3. 



, and Flux density for the frequencies 
6 and 6 cm respectively. 




02^920 OSESOa DS5S40 

RIGHT ASCENSJOM (J2Qaa) 



Fig. 1. 1.4 GHz radio image of A401 obtained from the NVSS (contours) overlaid on the ROSAT HRI X-ray image 
(colors) in the 0.1 — 2.4 keV band. The NVSS image has an angular resolution of 45'^ The first radio contour is drawn 
at 1.5 mJy/beam and the rest are spaced by a factor of ^/2. The high resolution observations at 6 cm (contours) 
overlaid on the DSS2 red plate (grey-scale) of the two radio galaxies are inset. A401A : Contour levels = 0.05 0.1 1 2 4 
6 mJy/beam. Angular resolution = 1.67^' x 1.42'^ (PA=— 45.2). Sensitivity = 0.015 mJy/beam. A401B : Contour levels 
= 0.05 0.1 0.2 0.5 1 mJy/beam. Angular resolution = 1.52'' x 1.38'' (PA=-19.8). Sensitivity = 0.016 mJy/beam. 



frequencies, calculated in three different locations of the 
source. 

3.2. Abell 2142 

The ROSAT HRI X-ray image presented in Fig. |4] shows 
the X-ray emission of the cooling flow cluster A2142. The 
cooling flow manifests itself with a cool X-ray peak as well 



as strong evidence for a centrally enhanced metal abun- 
dance (e.g. White et al. 1994, Peres et al. 1998). On the 
other hand optical (Oegerle et al. 1995) and X-ray (Henry 
& Briel 1996, Buote & Tsai 1996, Pierre & Starck 1998) 
features, such as the presence of an off-center, small scale 
structures and azimuthally asymmetric temperature vari- 
ations, are indications of a merger. Chandra data show 
that the gas temperature in the central region of A2142 
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8465 MHz 4535 MHz 

13 34 30 



^ 33 45 





02 58 32.5 32.0 31.5 31.0 30.5 30.0 29.5 29.0 28.5 28.0 
RIGHT ASCENSION (J2000) 



02 58 32.5 32.0 31.5 31.0 30.5 30.0 29.5 29.0 28.5 28.0 
RIGHT ASCENSION (J2000) 




5000 6000 7000 8000 
Frequency [MHz] 




5000 6000 7000 8000 
Frequency [MHz] 




5000 6000 7000 8000 
Frequency [MHz] 



Fig. 2. Source A401A. Top, left: Total intensity contours and polarization vectors at 3.6c7ti (8465 MHz). Contour 
levels are drawn at: -0.07 0.07 0.15 0.5 1 and 10 mJy/beam. Top, right: Total intensity contours and polarization 
vectors at 6cm (4535 MHz). Contour levels are drawn at: -0.1 0.1 0.5 1 and 10 mJy/beam. The angular resolution is 
3.0^' X 3.0'^ The lines give the orientation of the electric vector position angle (E-field) and are proportional in length 
to the fractional polarization {!" c^ 10%). Bottom: Trend of the fractional polarization as a function of the observing 
frequencies, at different source locations. The depolarization DP has been calculated as the ratio of the fractional 
polarization between the two extreme frequencies. 



(Markevitch et al. 2000) is non-isothermal. It appears that 
the central cooling flow has been disturbed but not de- 
stroyed by the merger. 

The presence of a halo in this cluster was suggested by 
Harris et al. (1977). Giovannini & Feretti (2000) confirmed 
the presence of diffuse emission, located around the bright- 
est cluster galaxy. This source, however, is much smaller 
than radio halos commonly found in merging clusters and 
it may be a mini-halo, like those found at the center of 
relaxed systems (see e.g. Gitti et al. 2007, Govoni et al. 
2009). 

In this work we have analyzed the radio galaxy 
labeled with A2142A in Fig. HI This narrow head 
tail radio galaxy has the head at a projected dis- 
tance of about 2.7' to the northwest of the clus- 
ter X-ray center and the tail elongated versus the 
west. The high resolution image at 6 cm shows 
the head, at the position R.A.(J2000)=15h58ml4.3s, 



Decl.(J2000)=27°16'19.5'', coincident with a galaxy. At 
6 cm the source is more than 40^' in size. In the source 
we do not distinguish the two jets. A possibility is that 
the source is narrow because of projection effects, i.e. the 
radio jets are seen edge-on. 

Fig. [5] shows the polarization images at 8465 MHz 
and 4535 MHz^ with an angular resolution of 3^' x 3" . At 
8465 MHz the polarization is significant (above Sa^poL) 
up to a distance of about 30'^ while at 4535 MHz the 
polarization signal appear below the noise level in some 
places along the tail. By considering all the pixels where 
the total intensity signal is above 5<ji at all frequencies, 
the mean fractional polarization is 16 ± 4% at 8465 and 
10 ±3% at 4535 MHz. The source is therefore depolarized 
< DP >= FPOU^s5MHz/FPOL8465MHz ^ 0.6. The de- 
polarization of the signal is also evident in the bottom 
panels of Fig. [5l where we show the trend of the fractional 
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Fig. 3. Source A401B. Top, left: Total intensity contours and polarization vectors at 3.6 cm (8465 MHz). Contour 
levels are drawn at: -0.07 0.07 0.15 0.5 1 and 3 mJy/beam. Top, right: Total intensity contours and polarization 
vectors at 6cm (4535 MHz). Contour levels are drawn at: -0.1 0.1 0.5 1 and 3 mJy/beam. The angular resolution is 
3.0'' X 3.0''. The lines give the orientation of the electric vector position angle (E-field) and are proportional in length 
to the fractional polarization (1" c^ 10%) Bottom: Trend of the fractional polarization as a function of the observing 
frequencies, at different source locations. The depolarization DP has been calculated as the ratio of the fractional 
polarization between the two extreme frequencies. 



polarization as a function of the observing frequencies, at 
different source locations. 



3.3. Abell 2065 

The ROSAT HRI X-ray image presented in Fig. [6] shows 
the X-ray emission in the 0.1 — 2.4 keV band of A2065. 
This is one of the galaxy clusters that make up the Corona 
Borealis super-cluster. The cluster X-ray surface bright- 
ness of A2065 is consistent with the optical galaxy surface- 
density profile, both showing elongation in the north- 
west - southeast direction. The cluster is a centrally con- 
densed system with two dominant galaxies separated on 
the plane of the sky by 17" and by :^ 600 km/ sec in red- 
shift (Postman et al. 1988). An ASCA X-ray observation 
of A2065 provides evidence for an ongoing merger, but 
also suggests two surface brightness peaks coincident with 
the two central galaxies, apparently survived the merger 
shock passage (Markevitch et al. 1999). A Chandra ob- 
servation (Chatzikos et al. 2006) reveals only one X-ray 



surface brightness peak, which is associated with the more 
luminous, southern galaxy. The gas related to that peak 
is cool and displaced slightly from the position of the 
galaxy. We note that the cluster temperature given in 
Table 1 (Markevitch 1998) is consistent with results from 
the ASCA and Chandra spectrum (White 2000, Ikebe et 
al. 2002, Chatzikos et al. 2006). On the other hand, David 
et al. (1993) derived a temperature of kT = S.A^H keV, 
from the Einstein MFC detector. 

No evidence for the presence of a cluster radio halo is 
available in the literature for A2065. 

In this work we have analyzed the radio galaxy la- 
beled with A2065A. It is located at a projected distance 
of about 13.6' from the cluster center. The high resolution 
image at 6 cm, overlaid on a R-band image taken with the 
Telescopio Nazionale Galileo (TNG; Fig. [6]), shows a radio 
galaxy of about 1' in size. It shows a FRII structure, with 
a core located at the position R.A.(J2000)=15h22m46.3s, 
Decl.(J2000)=27°55'40"and two opposite lobes, with evi- 
dence of hot-spots at the boundary of the lobes. 
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Fig. 4. 1.4 GHz radio image of A2142 obtained from the NVSS (contours) overlaid on the ROSAT HRI X-ray image 
(colors) in the 0.1 — 2.4 keV band. The NVSS image has an angular resolution of 45'^ The first radio contour is 
drawn at 1.5 mJy/beam and the rest are spaced by a factor of a/2- The high resolution observation at 6 cm (contours) 
overlaid on the DSS2 red plate (grey-scale) of the radio galaxy A2142A is inset. Contour levels = 0.05 0.1 0.2 0.5 1 2 
mJy/beam. Angular resolution = 1.39'^ x 1.18'' (PA=-57.5). Sensitivity = 0.014 mJy/beam. 



This radio galaxy does not show any optical counter- 
part in the DSS2 red image. The deeper, 30 min expo- 
sure R-band image we have taken with the TNG shows a 
faint source with a magnitude R=22.0±0.3. The small size 
of the optical counterpart in comparison with the clus- 
ter galaxies suggests that this object is a cluster back- 
ground source. Moreover, if we assume that the radio 
source is at the cluster redshift, the estimated absolute 
red magnitude (with no K-correction applied) would be 
Mr c^ 22.0 + 5 - 5 X IoqiqDl^pc) = -15.6. This is a 
very faint magnitude. Considering that the typical abso- 
lute magnitude of a FRII is Mr = —23.6 (e.g. Scarpa et 
al. 2002), the apparent magnitude of this radio source is 
explained if it is located at z > 1. 



Fig. shows the polarization images at 8465 MHz and 
4535 MHz^ with an angular resolution of 3'' x 3''. The po- 
larization has been detected in the lobes at each frequency. 
The two lobes have a similar fractional polarization. In the 
northern lobe the fractional polarization has a mean value 
of 16 ±5% both at 8465 MHz and 4535 MHz, while in the 
southern lobe the fractional polarization has a mean value 
of 19 ± 5% at 8465 MHz and 17 ± 5% at 4535 MHz. The 
absence of depolarization is also confirmed in the bottom 
panels of Fig. [71 where we show the trend of the fraction 
polarization as a function of the observing frequencies, 
calculated in three different locations of the source. 



3.4. Ophiuchus 

Ophiuchus is a nearby, hot cluster. Suzaku data (Fujita 
et al. 2008b) indicate that this cluster has a cool, dense 
core. We recently found at the cluster center a low surface 
brightness diffuse emission that it is classified has a mini- 
halo (Govoni et al. 2009, Murgia et al. 2009, Murgia et al. 
2010b). Moreover for this cluster there have been claims 
of detecting hard X-ray emission in excess to the ther- 
mal emission (Eckert et al. 2008, Nevalainen et al. 2009). 
Such emission may be of non-thermal origin, for example, 
Compton scattering of relativistic electrons by the cos- 
mic microwave background radiation (see e.g. Rephaeli et 
al. 2008, Petrosian et al. 2008, and references therein for 
recent reviews), but alternative explanations have been 
put forward (Profumo 2008, Perez-Torres et al. 2009, 
Colafrancesco & Marchegiani 2009). 

Fig.[8]shows the ROSAT HRI X-ray image of the clus- 
ter overlaid on the NVSS image. In this work we have ana- 
lyzed the radio galaxies labeled with OPHIB and OPHIC. 
In the original sets of data we observed also OPHIA (see 
Fig [8]) , but this source is fully resolved out at both wave- 
lengths. The lack of detection indicates that the OPHIA 
source is too extended to be imaged at a such small angu- 
lar scale. 

The source OPHIB is located to the northwest 
of the cluster center at a projected distance of 
about 14.5'. The high resolution image at 6 cm over- 
laid on the DSS2 red plate shows a source of 
about 45" in size with a NAT head-tail morphology. 
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Fig. 5. Source A2142A. Top, left: Total intensity contours and polarization vectors at 3.6 cm (8465 MHz). Contour 
levels are drawn at: -0.06 0.06 0.15 0.5 and 3 mJy/beam. Top, right: Total intensity contours and polarization vectors at 
6 cm (4535 yiHz). Contour levels are drawn at: -0.1 0.1 0.5 1 and 3 mJy/beam. The angular resolution is 3.0''x3.0'^ The 
lines give the orientation of the electric vector position angle (E-field) and are proportional in length to the fractional 
polarization (1^^ c^ 10%). Bottom: Trend of the fractional polarization as a function of the observing frequencies, 
at different source locations. The depolarization DP has been calculated as the ratio of the fractional polarization 
between the two extreme frequencies. 



The core at the position R.A.(J2000)=17hllm55.4s, 
Decl.(J2000)=-23°09'42.5'' is coincident with a bright 
galaxy, while the tail is oriented away from the core in 
the northwest direction. 

The NAT head-tail OPHIC is located at a projected 
distance of about 7.6' to the southwest of the cluster center 
with the head pointed versus the cluster center and the tail 
elongated versus south. The high resolution image at 6 cm 
shows the core, at the position R.A.(J2000)=17hl2m09.0s, 
Decl.(J2000)=-23°28'26.5'', coincident with a galaxy. At 
6 cm the source is about 40^' in size. 

Figs. [9] and [10] show the polarization images at 8465 
yiHz and 4835 Mi^z, with an angular resolution of 3^' x 3'^ 

OPHIB is polarized at both frequencies up to a dis- 
tance of about 45'^ By considering all the pixels where 
the total intensity signal is above 5<ji at all frequencies, 
the source has a mean fractional polarization of 18 ± 5% 
at 8465 and 22 ± 5% at 4835 MHz. Within the errors the 
two fractional polarizations are compatible therefore in 
OPHIB we do not detect significant depolarization. In the 
bottom panels of Fig.[9]we show the trend of the fractional 
polarization as a function of the observing frequencies, at 
different source locations. 

The fractional polarization in different locations of the 
source OPHIC is indicated in Fig. [JOl At 8465 MHz the 



polarization is detected only in a few places and the error 
on the fractional polarization is rather high. In most of 
the source the fractional polarization at 4835 is below the 
3crFPOL level. The only patch where a significant polar- 
ization is detected has a fraction polarization compatible, 
within the errors, with that at 8465 MHz. 

4. Rotation measure images 

The presence of a magnetized and ionized screen between 
the observer and a radio source changes the properties of 
the incoming polarized emission. As it passes through the 
magnetized plasma, the polarized synchrotron radiation 
undergoes the following rotation of the plane of polariza- 
tion: 



^OhsW = ^\nt^>? xRM 



(1) 



where ^o6s(A) is the observed polarization angle at a 
wavelength A and "^int is the intrinsic polarization angle. 
The rotation measure RM is related to the plasma 
thermal electron density, ng, and magnetic field along the 
line-of-sight, ^m, by the equation: 



RM = 812 / UeBudl lad m~ 



(2) 
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Fig. 6. 1.4 GHz radio image of A2065 obtained from the NVSS (contours) overlaid on the ROSAT HRI X-ray image 
(colors) in the 0.1 — 2.4 keV band. The NVSS image has an angular resolution of 45'^ The first radio contour is drawn 
at 1.5 mJy/beam and the rest are spaced by a factor of a/2- The high resolution observation at 6 cm (contours) overlaid 
on a TNG R-band image (grey-scale) of the radio galaxy A2065A is inset. Contour levels = 0.05 0.08 0.12 0.15 0.2 0.5 
12 3 5 mJy/beam. Angular resolution = 1.24'' x 1.18'' (PA=-45.5). Sensitivity = 0.014 mJy/beam. 

Table 4. Rotation measure results. 



Source 


Distance 


N 


<Errfit> 


(RM) 


ctrm 




(kpc) 




(rad/m^) 


(rad/m^) 


(rad/m^) 


A401A 


440 


34 


30 


222±20 


113±14 


A401B 


730 


23 


41 


107=bl5 


74zbl4 


A2142A 


270 


8 


36 


-461il=82 


230±63 


A2065A 


1120 


12 


29 


-9±16 


48=bl5 


OPHIB 


480 


15 


44 


-64±22 


74±19 



Col. 1: Source; Col. 2: Projected distance from the X-ray center; 
Col. 3: Number of beams over which RMs are computed; 
Col. 4: Mean value of the RM fit error; 
Col. 5: Mean of the RM distribution; 
Col. 6: RMS of the RM distribution; 



where B\\ is measured in jaG, 
depth of the screen in kpc. 



Up in cm and L is the 



Following the definition, the RM images were obtained 
by performing a fit of the 4 polarization angle images, at 
each pixel, as a function of A^ (see Eq(T]), by using the 
software FARADAY by Murgia et al. (2004). Given as 
input the images of Q and [/, at each frequency, the soft- 
ware produces the RM and the intrinsic polarization angle 
^int images, both with relative error images. To improve 
the RM image, the software can be iterated in several 
self-calibration cycles. In the first cycle only pixels with 
the highest signal-to-noise ratio are fitted. In the next cy- 
cles the algorithm uses the RM information in these high 
signal-to-noise pixels to solve the UTr-ambiguity in adjacent 



pixels of lower signal-to- noise, in a similar method used in 
the PACERMAN algorithm by Dolag et al. (2005b). 

We derived the RM images, at 3^' resolution, of the 
five polarized sources. They were calculated only in those 
pixels in which the following three conditions were sat- 
isfied: the total intensity signal at 8465 MHz was above 
3(7/ , the error in the polarization angle at each frequency 
was lower than 15° , and the resulting RM error was lower 
than 60 rad/m^. We checked the agreement between the 
final RM images obtained with FARADAY and those ob- 
tained with the algorithm PACERMAN by Dolag et al. 
(2005b). 

Figs. 11 — 15 show the resulting RM images for the 
sources A401A, A401B, A2142A, A2065A and OPHIB, re- 
spectively. For each source the total intensity contours at 
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Fig. 7. Source A2065A. Top, left: Total intensity contours and polarization vectors at 3.6c7ti (8465 yiHz). Contour 
levels are drawn at: -0.06 0.06 0.1 0.5 and 3 mJy/beam. Top, right: Total intensity contours and polarization vectors 
at 6cm (4535 MHz). Contour levels are drawn at: -0.1 0.1 0.15 3 and 6 mJy/beam. The angular resolution is 3.0'^ x 
3.0'^ The lines give the orientation of the electric vector position angle (E-field) and are proportional in length to 
the fractional polarization (1" c^ 10%). Bottom: Trend of the fractional polarization as a function of the observing 
wavelengths, at different source locations. The depolarization DP has been calculated as the ratio of the fractional 
polarization between the two extreme frequencies. 



8465 MHz are overlaid on the RM images. We can char- 
acterize the RM distribution in terms of a mean {{RM)) 
and root mean square {(Jrm)- The histograms of the RM 
distribution are also shown. In addition, a few plots show- 
ing the position angle ^o6s as a function of A^ at different 
source locations are presented. 

In Table 4, for each radio galaxy we report its pro- 
jected distance from the cluster X-ray center, the num- 
ber of beams over which RMs are computed, the mean 
value of the RM fit error, the mean {{RM)) and root 
mean square {(Trm) of the RM distribution. The errors 
on the {RM) and cfrm quantities have been computed 
with Montecarlo simulations, considering both the uncer- 
tainty related to the presence of the statistical error and 
the fit erroo. 



^ Indeed, our data suffer from two kinds of uncertainty. The 
first source of error is related to the fit of the A^-law and takes 
into account the presence of errors in the measurements of Q 



The values of (RM) given in Table 4 are not cor- 
rected for the contribution of the Galaxy. We determined 
the Galactic contribution for each cluster of our sample 
according to the values reported in the RM catalog by 
Taylor et al. (2009). We used the distance- weighted mean 
of the RMs of all sources located within a radius of 3 deg. 
The RM Galactic contribution in the region occupied by 
our cluster sample has been found to be rather small. 
About -6 rad/m^ for A401, +16 rad/m^ for A2142, and 



and U . The fit error has the effect to broad the histogram of 
the RM distribution and thus to increase the observed grm- 
While the second uncertainty is the statistical error. The sta- 
tistical errors for {RM) and for grm is given by (trm /^/N and 
ctrm/V'^N respectively, where N is the number of beams over 
which the RM has been computed. Through the Montecarlo 
procedure we checked that the fit error is negligible with re- 
spect to the statistical error. 
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Fig. 8. 1.4 GHz radio image of Ophiuchus obtained from the NVSS (contours) overlaid on the ROSAT HRI X-ray 
image (colors) in the 0.1 — 2.4 keV band. The NVSS image has an angular resolution of 45'^ The first radio contour is 
drawn at 1.5 mJy/beam and the rest are spaced by a factor of a/2. The high resolution observations at 6 cm (contours) 
overlaid on the DSS2 red plate (grey-scale) of the two radio galaxies are inset. Contour levels = 0.1 0.15 0.3 0.5 0.8 
1 mJy/beam. OPHIB : Angular resolution = 1.09'' x 0.75'' (PA=-87.0). Sensitivity = 0.038 mJy/beam. OPHIC : 
Angular resolution = 1.09" x 0.76" (PA=-87.7). Sensitivity = 0.030 mJy/beam. 



+10 rad/m^ for A2065. While the Galactic RM is 
rad/m^ for Ophiuchus. 



-29 



In order to verify the polarization angle linearity with 
A^, we selected some pixels in different locations of the 
sources. The data are generally quite well represented by 
a linear A^ relation. The polarization angle trend with A^ 
and the absence of depolarization (at least in the case 
of A401, A2065, and Ophiuchus) would argue in favor of 
an RM produced by an external Faraday screen (Burn 
1966, Laing 1984, Laing et al. 2008). However, some au- 
thors have suggested the possibility that the RM observed 
in radio galaxies is not associated with the foreground 
intracluster medium, but may arise locally to the radio 
source (Bickneh et al. 1990, Rudnick & Blundeh 2003). 
A much larger A^ range would be needed to unambigu- 
ously exclude the presence of an internal Faraday screen. 
However, although in this work we cannot address the ori- 
gin of the screen properly (internal, external, or both), 2D 
Montecarlo simulations presented in the Appendix show 
that both RM and polarization data can be explained by 
the presence of a foreground Faraday screen only. In agree- 
ment with the 2D Montecarlo simulations, in the following 
we assume that the internal contribution is negligible and 
that the major factor responsible for the Faraday rotation 
is given by the external intracluster medium (e.g. Govoni 
& Feretti and references therein). In this case, the struc- 
ture of the RM images and the values of the cfrm and 
(RM) can give information on the intracluster magnetic 



field structure. In particular, the RM dispersion of the his- 
tograms and the presence of small scales structure in the 
RM images can be explained by the fact that the cluster 
magnetic field fluctuates on scales smaller than the size of 
the sources. On the other hand, the cases of RM distribu- 
tions (once corrected for the Galactic contribution) with a 
non-zero mean indicate that the magnetic field fluctuates 
also on scales larger than the radio sources. Therefore, in 
order to interpret correctly the RM data, magnetic field 
fluctuations in a wide range of spatial scales should be 
considered. 

The radio galaxies analyzed in this work are located at 
a projected distance from the cluster center which range 
from 270 kpc (A2142A) up to 1120 kpc (A2065A). 

The case of A401 is very interesting because it allows 
the RM to be sampled along two different lines-of-sight 
within the same cluster. As expected in the case of the 
RM being due to an external Faraday rotation due to 
the cluster magneto-ionized medium, the innermost source 
(A401A) has a higher ctrm and |(RM)| than the external 
one (A401B). The RM images reveal patchy structures 
with RM fluctuations down to scales of a few kpc. 

These small scale RM structures are well visible also in 
the other radio galaxies analyzed here. In A401A, A401B 
and A2142 large patches of fairly constant RM are also 
visible, indicating a magnetic structure fluctuating on 
large scales. We note however that the sampling of large 
scale fluctuations is rather poor. 
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Fig. 9. Source OPHIB. Top, left: Total intensity contours and polarization vectors at 3.6 cm (8465 MHz). Contour 
levels are drawn at: -0.15 0.15 0.5 1 2 and 3 mJy/beam. Top, right: Total intensity contours and polarization vectors 
at 6cm (4835 MHz). Contour levels are drawn at: -0.35 0.35 0.7 1 2 and 3 mJy/beam. The angular resolution is 
3.0'' X 3.0''. The lines give the orientation of the electric vector position angle (E-field) and are proportional in length 
to the fractional polarization (1" ~ 10%). Bottom: Trend of the fractional polarization as a function of the observing 
wavelengths, at different source locations. The depolarization DP has been calculated as the ratio of the fractional 
polarization between the two extreme frequencies. 



We note that the location of Ophiuchus near the 
Galactic center makes this target not optimal for a 
magnetic field structure investigation through the RM. 
Therefore, the RM image of OPHIB should be considered 
with caution. 

A2065A is a nice example of a high resolution RM 
image in an extended background source. Despite its dis- 
tance from the cluster center, its importance is due to fact 
that, at least in this case, we can exclude that the patchy 
structures seen in the RM image are due to the turbu- 
lence locally created by the motion of the galaxy in the 
intracluster medium. 



5. Magnetic field - gas temperature connection 

The new RM data presented here increase the statistic 
of the RM for radio sources in hot clusters. Therefore, 
by comparing these new data with RM information taken 
from the literature it is now possible to investigate the 
connection between the magnetic field strength and the 



temperature of the intracluster medium. For a detailed 
statistical analysis, we formed a sample of 12 galaxy clus- 
ters for which high-quality RM data are available. 

The cluster X-ray properties of the sample are pre- 
sented in Table 5. For the distribution of the thermal elec- 
tron gas density we assumed the standard /3-model profile 
(Cavaliere & Fusco-Femiano 1976): 



ne(r)=no(l + rVr,')-'^/' 



(3) 



where no and Tc are the central electron density, and the 
cluster core radius, respectively. The X-ray parameters re- 
ported in Table 5 are taken from the literature and cor- 
rected for the cosmology adopted in this paper. The mean 
cluster density calculated by integrating the /3-model pro- 
file over a sphere of IMpc in radius is also indicated. 

In the sample the temperature of the clusters is in the 
range T c:^ 3 — 10 keV, while the mean /3- model parameters 
are: < (3 >=0.7, < Tc >=285 kpc, and < no >= 4.4 x 10"^ 
cm~^. In the cluster sample there are a total of 29 radio 
galaxies, for which high-quality RM data are available. 
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Fig. 10. Source OPHIC. Left: Total intensity contours and polarization vectors at 3.6 c?n (8465 MHz). Contour levels 
are drawn at: -0.1 0.1 0.5 1 and 2 mJy/beam. Right: Total intensity contours and polarization vectors at 6 cm (4835 
MHz). Contour levels are drawn at: -0.3 0.3 0.7 2 and 3 mJy/beam. The angular resolution is 3.0'^ x 3.0'^ The lines 
give the orientation of the electric vector position angle (E-field) and are proportional in length to the fractional 
polarization (1'' c^ 17%). 



For each radio source, in Table 6 we report the projected 
distance from the cluster X-ray center, and the ctrm taken 
from the literature. We included in this analysis only those 
clusters hosting radio sources for which high-resolution, 
extended RM images are available. We did not consider 
RM images of sources located at the center of cooling 
core clusters (e.g. Taylor et al. 2002). Indeed in our list 
there are a few cooling flows clusters (A2142, Ophiuchus, 
A2382), but in these cases the analyzed sources are lo- 
cated at a large projected distances from the cluster cen- 
ter, therefore the contribution of the dense cool core to 
the RM is likely to be negligible. 

In Fig. [16] (top) we analyze the trend of ctrm^ redshift 
corrected by multiplying by a factor of (1 + z)^, as a func- 
tion of distance from the cluster X-ray center. The differ- 
ent symbols represent different cluster temperature. It is 
evident a clear broadening of the aRM toward small pro- 
jected distances, consistently with an excess of Faraday ro- 
tation due to the magnetized intracluster medium. In ad- 
dition, a positive trend is found between the cluster tem- 
perature and the aRM- In particular, for a fixed projected 
distance from the cluster center, clusters with higher tem- 
perature show a higher ctrm- 

The interpretation of ctrm in a cluster depends on sev- 
eral assumptions, including the model for the X-ray emit- 
ting gas and the magnetic field structure. To try to explain 
the trend between the cluster temperature and the aRM 
seen in Fig. [16] (top), for simplicity, we can consider an an- 
alytical formulation based on the approximation that the 
magnetic field is constant along the cluster and tangled 



on a single scale A^. In this simple ideal case, the screen 
is made of cells of uniform size, electron density and mag- 
netic field strength, but with a field orientation at random 
angles in each cell. The observed RM along any given line 
of sight is then generated by a random walk process involv- 
ing a large number of cells of size A^. The distribution of 
the RM is Gaussian with zero mean, and a variance given 
by: 



1^ 



(4) 



In this formulation, by considering a density distribu- 
tion which follows a /3-profile, the following relation (e.g. 
Lawler & Dennison 1982, Tribble 1991, Feretti et al. 1995, 
Felten 1996) for the RM dispersion is obtained by inte- 
grating Eq. m 



CTRMir) 



(1^^2/^2)(6/3-l)/4' 



/ r(3/3-0.5) 

r(3/3) 



(5) 



where B = V3B\\ and F is the Gamma function. The 
constant K depends on the integration path over the gas 
density distribution: K = 624, if the source lies completely 
beyond the cluster, and i^ = 441 if the source is halfway 
the cluster. Therefore, aRM depends on the possible fluc- 
tuations in no, 5, v^c and a/Ab- Since the linear depen- 
dence on no and B strongly affect aRM, in the following 
we concentrate on these two parameters. For most clusters 
of our sample, measurements of the magnetic field scale 
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Fig. 11. Rotation measure image of the radio galaxy A401A. The angular resolution is 3.C x 3.0'^ The contours 
refer to the total intensity image at 8465 M^Hz. Contour levels are drawn at: 0.07 0.15 0.5 1 and 5 mJy/beam. The 
histogram shows the rotation measure distribution for all significant pixels. The plots show the position angle ^ohs as 
a function of A^ at different source locations. 



are not available from the literature. However, a discus- The trend seen in the top of Fig. [161 can be explained 

sion regarding the effect of A^ on grm is reported in the by a gas temperature — gas density relation, or by a gas 
Appendix. temperature — magnetic field strength relation. 
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Fig. 12. Rotation measure image of the radio galaxy A401B. The angular resolution is 3.0^' x 3.0'^ The contours 
refer to the total intensity image at 8465 M.liz. Contour levels are drawn at: 0.07 0.15 0.5 1 and 3 mJy/beam. The 
histogram shows the rotation measure distribution for all significant pixels. The plots show the position angle ^ohs as 
a function of A^ at different source locations. 



In the bottom panels of Fig.[T6lwe show the predictions sample (/3=0.7, Vc = 285 kpc). On the left bottom panel, 

of the analytical formulation (Eq.[5]) obtained by assuming we fixed the central gas density to no = 4.4 x 10~^ cvrT^ 

a magnetic field scale A^ = 10 kpc, and by keeping fixed (the mean value of the cluster sample) and we present the 

the X-ray parameters to the mean values of the cluster expectations for three different magnetic field strengths 
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Fig. 13. Rotation measure image of the radio galaxy A2142A. The angular resolution is 3.0^' x 3.0'^ The contours 
refer to the total intensity image at 8465 MHz. Contour levels are drawn at: 0.06 0.15 0.5 1 and 3 mJy/beam. The 
histogram shows the rotation measure distribution for all significant pixels. The plots show the position angle ^o6s as 
a function of A^ at different source locations. 



{B = 1,2,4 fiG). On the right bottom panel, we fixed produced quite well. Therefore, the higher aRM in hotter 
B = 3 jiiG and we present the expectations for three dif- clusters may be explained if they have a higher magnetic 



ferent central gas densities (no = 1.5,3,6 x 10 



). field strength and/or if they have a higher gas density. 



In both cases the statistical trend of the data can be re- 
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Fig. 14. Rotation measure image of the radio galaxy A2065A. The angular resolution is 3.C x 3.0'^ The contours 
refer to the total intensity image at 8465 M.liz. Contour levels are drawn at: 0.06 0.1 0.5 and 3 mJy/beam. The 
histogram shows the rotation measure distribution for all significant pixels. The plots show the position angle ^o6s as 
a function of A^ at different source locations. 



It is not easy to disentangle the two effects; however, the temperature of the clusters. Clusters with a higher 

in the following we show that hotter clusters have higher temperature show a higher gas density. The difference of 

(yRM mostly because of their higher gas density. In Fig. about a factor of four in the mean density between cooler 

[T71 (left), we show the mean gas density plotted versus and hotter clusters may explain the difference of about a 
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Fig. 15. Rotation measure image of the radio galaxy OPHIB. The angular resolution is 3.C x 3.0'^ The contours refer 
to the total intensity image at 8465 Mi^z. Contour levels are drawn at: 0.15 0.5 1 2 and 3 mJy/beam. The histogram 
shows the rotation measure distribution for all significant pixels. The plots show the position angle ^ohs as a function 
of A^ at different source locations. 



factor of four in grm between cooler and hotter clusters. In Fig. [TTl (right), we show the core radius plotted ver- 

Therefore, this analysis does not confirm a strict link be- sus the temperature of the clusters. No obvious trend is 

tween the magnetic field strength and the gas temperature present between the temperature and core radius of the 

of the intracluster medium. clusters. 
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Fig. 16. Top: Dispersion of the rotation measure distribution as a function of the projected distance from the cluster 
X-ray center. The different symbols represent the cluster temperature taken from the literature (red > 7 keV, green 4 — 7 
keV, blue < 4 keV). Bottom: Prediction of the analytical formulation by assuming a magnetic field scale A^ = 10 kpc, 
and by keeping fixed the X-ray parameters to the mean values of the cluster sample /3 = 0.7 re = 285 kpc. On the 
left no = 4.4 x 10~^ cm~^ has been fixed and we present the expectations for three different magnetic field strengths 
{B = 1,2,4 /iG). On the right B = 3 fiG has been fixed and we present the expectations for three different central 



gas densities (no = 1.5,3,6 x 10 



')■ 



An alternative method to investigate a possible con- 
nection between the magnetic field strength and the gas 
temperature of the intracluster medium is the analysis 
of the CF^M — Sx correlation. The cluster X-ray surface 
brightness Sx is given by: 



Sx oc 



]^/Tdl 



LOS 



(6) 



The X-ray surface brightness (see Table 6) has been 
calculated by analyzing pointed ROSAT PSPC observa- 
tions in the 0.1 - 2.4 keV band. ROSAT PSPC pointed 
observations were not available in the case of 3C129 
and A2065 therefore we used the images taken from the 
ROSAT All Sky Survey, in the same band. The X-ray sur- 
face brightness has not been calculated in the location 
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of the sources 5C4.127 and 5C4.42 (in the Coma clus- 
ter) because they are located outside the field of view of 
the ROSAT detector. The X-ray surface brightness (in 
Cts s~^pixel~^) has been derived by averaging the X-ray 
count rates in an annulus comprising the position of the 
radio galaxies. We then corrected the count rates for the 
background. The background has been measured well out- 
side the X-ray-bright cluster region. X-ray point sources 
have been masked out in this analysis. We finally con- 
verted the X-ray surface brightness from Cts s~^ pixel" ^ 
to erg s~^cm~^sterad~^ by using the PIMMS softwarqj. In 
this conversion the X-ray emission was approximated by 
a Raymond- Smith model with the mean cluster tempera- 
ture and photoelectric absorption column density given in 
the literature (see Table [5]). 

In Fig. [18] we report the aRM — Sx plot showing the 
trend of (Trm values as a function of the cluster X-ray sur- 
face brightness Sx calculated in the corresponding source 
location. In the plot, Sx has been corrected for the de- 
pendence on the temperature by multiplying by a factor 
of T~^-^. In addition, Sx has been corrected for the cos- 
mological dimming of the surface brightness by dividing 
by a factor of (1 -\- z)^. The new RM data presented here 
(A401A, A401B, OPHIB, A2142A, A2065A) confirm the 
o'RM — Sx relation previously known in the literature (e.g. 
Dolag et al. 2001, Weratschnig et al. 2008). The differ- 
ent symbols represent different cluster temperature. As 
predicted by cosmological magneto-hydrodynamic simu- 
lations (Dolag et al. 1999), galaxy clusters may have dif- 
ferent central magnetic field strengths depending on their 
dynamical state, leading to an offset of the aRM — Sx 
relation, along the aRM axis, depending on the cluster 
temperature. An offset of the aRM — Sx relation, along 
the aRM axis, depending on the cluster temperature if 
present is comparable to the intrinsic scatter of the corre- 
lation. Therefore, a possible connection between the mag- 
netic field strength and the gas temperature, if present, 
is very weak. In addition, although the statistics is still 
limited, we note that in the aRM — Sx relation clusters 
with extended diffuse emission (e.g. Coma, A2255, A401) 
seem to show a behavior similar to that of clusters without 
diffuse synchrotron emission. 

In conclusion, present data does not permit to establish 
a strict link between the magnetic field strength and the 
gas temperature of the intracluster medium. 

6. Conclusions 

We performed Very Large Array observations at 3.6 cm 
and 6 cm of two radio galaxies located in A401 and 
Ophiuchus, a radio galaxy in A2142, and a radio galaxy 
located in the background of A2065. For most of them we 
obtained detailed, high resolution (:^ 3^'), rotation mea- 
sure images. These RM images reveal patchy structures 
with RM fluctuations down to scales of a few kpc. Under 
the assumption that the radio galaxies themselves have no 



effect on the measured RMs, these structures indicate that 
the intracluster magnetic fields fluctuate down to such 
small scales. 

These new RM data are compared with RM infor- 
mation taken from the literature for which similar high- 
resolution, extended RM images are available. 

A positive trend is found between the cluster temper- 
ature and the dispersion of the RM distributions (aRM)- 
The correlation found between the gas density and the 
gas temperature indicates that in our sample hotter clus- 
ters have higher aRM mostly because of their higher gas 
density. Moreover, although the previously known relation 
between the clusters X-ray surface brightness (Sx) at the 
radio galaxy location and aRM is confirmed, a possible off- 
set of the aRM — Sx relation along the aRM axis, depend- 
ing on the cluster temperature, if present is comparable 
to the intrinsic scatter of the correlation. 

Therefore, in view of the current data it is not possi- 
ble to establish a strict link between the magnetic field 
strength and the gas temperature of the intracluster 
medium. 
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Appendix A: Rotation measure and depolarization 
in a foreground Faraday screen 

Under the assumption that the RM spatial variations 
seen in the presented work are fully explained by external 
patchy screens, we investigated if these RM structural in- 
formation are consistent with the observed depolarization 
levels at longer wavelengths. Indeed, the most direct con- 
sequence arising from RM structures which fluctuate on 
small scales is the beam depolarization effect. 

We performed this analysis through the software 
FARADAY (Murgia et al. 2004). In particular we per- 
formed 2D Montecarlo simulations by assuming for the 
rotation measure a simple power law power spectrum of 
the type: 



\RMk 



(A.1) 



http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html 



For each source we simulated RM images by keeping 
fixed n = 11/3, corresponding to the Kolmogorov power 
law index for the parent magnetic field power spectrum, 
and by varying the minimum and maximum scale of the 
RM fluctuations (A^^^, ^max)- The computational grid 
was 1024 X 1024 pixel^ and the pixel-size was fixed to 1 
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Fig. 17. Left: Plot of the mean gas density versus the temperature of the clusters in the sample. The mean cluster 
density has been calculated by integrating the /3-model profile over a sphere of IMpc in radius. Right: Plot of the core 
radius versus the temperature of the clusters in the sample. The different symbols represent the cluster temperature 
taken from the literature (red > 7 keV, green 4 — 7 keV, blue < 4 keV). 

Table 5. X-ray data of the cluster sample. 



Cluster 


z 


kpcr 


T 

[keV] 


nH 


P 


kpc 


no 
10-^ cm-3 


<ne> 

10"^ cm-2 


Reference 


A514 


0.0714 


1.34 


3.8 


3.64x10^^ 


0.6 


417 


0.5 


0.15 


1,7 


Coma 


0.0232 


0.46 


8.38 


8.54x10^^ 


0.654 


245 


3.5 


0.44 


2,8 


A2255 


0.0806 


1.50 


6.87 


2.49x10^° 


0.797 


438 


2.1 


0.48 


3,8 


A400 


0.0240 


0.48 


2.31 


8.33x10^0 


0.534 


110 


2.4 


0.14 


2,8 


A2634 


0.0312 


0.61 


3.7 


4.79 xlO^^ 


0.640 


261 


2.8 


0.40 


2,8 


A119 


0.0441 


0.86 


5.8 


3.51x10^0 


0.675 


362 


1.8 


0.39 


4,8 


3C129 


0.0223 


0.44 


5.6 


5.96x10^1 


0.601 


226 


2.1 


0.27 


5,8 


A2382 


0.0618 


1.18 


2.9 


3.99x10^° 


0.9 


373 


1.2 


0.17 


6,9 


A401 


0.074 


1.39 


8.3 


9.88x10^° 


0.613 


177 


7.0 


0.60 


4,8 


A2142 


0.091 


1.67 


8.8 


3.78x10^0 


0.591 


114 


18.7 


0.86 


4,8 


A2065 


0.073 


1.37 


5.4 


3.04x10^° 


1.162 


507 


2.3 


0.39 


4,8 


Ophiuchus 


0.028 


0.55 


10.26 


1.93x10^^ 


0.747 


199 


8.0 


0.53 


2,8 



Col. 1: Cluster name; Col. 2: Redshift; Col. 3: Angular to linear conversion; Col. 4: Cluster temperature; 

Col. 5: Galactic absorption taken from the Leiden/Argentine/Bonn (LAB) Survey of Galactic HI (Kalberla et al. 2005); 

Col. 6, 7: /3-model parameters {(3 and core radius); Col. 8: Central gas density; 

Col. 9: Mean gas density calculated by integrating the ^-model profile over a sphere of IMpc in radius; 

Col. 10: Temperature and /3-model references: 

1 Weratschnig et al. (2008), 2 Fukazawa et al. (1998), 3 White (2000), 4 Markevitch (1998), 

5 Edge & Stewart (1991), 6 Ebeling et al. (1996), 

7 Govoni et al. (2001), 8 Chen et al. 2007, 9 Guidetti et al. (2008) 



kpc. The resulting field of view was then ^ 1024x1024 
kpc^, that was enough to recover the projected size of 
the sources and to properly sample the RM power spec- 
trum scales. We then added to the simulated RM images a 
Gaussian noise having dnoise = Ervfu^ in order to mimic 
the effect of the noise in the observed RM images and 
we then added an RM offset corresponding to the RM 
Galactic contribution in the sources location. Finally, the 
simulated RM images were convolved with a Gaussian 
function having FWHM equal to the beam of the obser- 
vations, and blanked as the observations. 



We performed several 2D simulations of synthetic RM 
images to evaluate the range of Kmin--^max which could 
explain the observed RM images. Fig. A.l (left panels) 
shows some examples of simulated RM images, which we 
found to reproduce quite well the observed RM distribu- 
tions. 

To quantify the two dimensional fluctuations of RM 
and to compare the simulations with the observations, we 
used the structure function as a statistical tool. The RM 
structure function is defined by 

S{dx,dy) = {[RM{x,y)-RM{x^dx,y^dy)f)^x^y),{k.2) 
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Fig. 18. Dispersion of the rotation measure distribution as a function of the X-ray surface brightness of the intracluster 
gas in the source location. The different symbols represent the cluster temperature taken from the literature (red > 7 
keV, green 4-7 keV, blue < 4 keV). 



where = (){x,y) indicates that the average is taken over all 
the positions {x^y) in the RM image. Blank pixels were 
not considered in the statistics. The structure function 
Sir) was then computed by radially averaging S{dx^ dy) 
over regions of increasing size of radius r = ^J dx^ + dy^ . 
In the middle panels of Figs. A.l we compare the observed 
(black points) and the simulated (grey-dashed) RM struc- 
ture function. 

In the case of an external Faraday screen, a reduction 
of the fractional polarization at longer wavelengths is ex- 
pected if the minimum scale of the RM fluctuation (A^^^) 
extend to small scales (beam depolarization). Once a good 
representation of the RM structures have been found, we 
investigated the expected beam depolarization to check 
the consistency with the observed polarization at different 
wavelengths. In particular, we produced, at each observ- 
ing frequency, the expected Q and U images correspond- 
ing to the simulated RM^ and we taken into account the 
beam depolarization effects by convolving Q and \J to the 
appropriate resolution. The predicted and observed mean 
degrees of polarization are then compared. 

For each source we found a combination of model pa- 
rameters ts^min^ ^max that givcs a reasonable representa- 
tion of the RM structure across the source as well as of the 



depolarization at longer wavelengths. The values of l^rain 
and ts^max we found for these sources are consistent with 
those found in the literature for a few other galaxy clus- 
ters (e.g. Murgia et al. 2004, Govoni et al. 2006, Guidetti 
et al. 2008, Bonafede et al. 2010). Therefore, even if we 
cannot exclude a possible presence of an internal contri- 
bution to the observed RM images, the 2D Montecarlo 
simulations indicate that both RM and polarization data 
can be explained by the presence of a foreground Faraday 
screen. 

The different magnetic field scales A^ of the clusters 
can be compared by calculating the magnetic field auto- 
correlation length, which takes into account the minimun 
and maximum scale of fluctuations (A^^^, ts.j^ax\ and the 
index of the power spectrum (n). Following the magnetic 
field autocorrelation length formula presented by EnBlin 
& Vogt (2003), for the new data analyzed here we found 
a magnetic field autocorrelation length variation of a fac- 
tor of 9. This difference may influence the grm values of 
about factor of 3 (grm oc \/Kb)- However, we note that 
the exact determination of the magnetic field power spec- 
trum is beyond the scope of this paper, and the numbers 
indicated in the caption of Fig. A.l should be considered 
as indicative. 
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Fig. A.l. Left: Examples of simulated RM images for the sources A401A, A401B, A2142, A2065, OPHIB. The 
simulated images have the same colour scale and resolution of the observed RM images. They have been obtained by 
simulating a RM power spectrum with a spectral index n=ll/3 and a minimum scale of the RM fluctuations b^rnin—^ 
kpc. The maximum scale of the RM fluctuations is A^ax=205 kpc for A401A and k2\^2^ A^aaj=128 kpc for A401B, 
Amax=15 kpc for OPHIB and A2065. Middle: Structure functions of the simulated (grey) and observed (black points) 
RM images. Right: Depolarization plots of the simulations (grey) and observations (black points). The depolarization 
is defined as DP{y) = FPOL^/FPOLsag^uhz' 
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Table 6. RM and X-ray data. 



Cluster 


Source 


Distance 


CTRM 


Sx [0.1 - 2.4] keV 


Reference 






[kpc] 


[rad/m2] 


10"'^ [erg/ s cm^ sterad] 




A514 


A514B2 


232 


63 


5.75zb0.13 






A514D 


509 


47 


1.56±0.03 






A514A 


1556 


48 


<1.79 






A514E 


1654 


47 


<1.79 






A514C 


1822 


24 


<1.79 




Coma 


5C4.85 


51 


303 


86.74±0.48 


2 




5C4.81,NGC4869 


124 


166 


65.22±0.26 


2 




5C4.74 


372 


154 


16.37±0.08 


2 




5C4.114 


532 


16 


6.76±0.05 


2 




5C4.127 


919 


65 


0.70±0.01 


2 




5C4.42 


1250 


56 


- 


2 




5C4.152 


1489 


37 


- 


2 


A2255 


1712.4+6401 


279 


79 


15.78±0.22 


3 




J1713.5+6402 


444 


59 


7.86±0.11 


3 




J1713.3+6347 


1497 


42 


0.33zb0.01 


3 


A400 


3C75 





100 


8.69±0.15 


4 


A2634 


3C465 





120 


10.77±0.30 


4 


A119 


0053-015 


124 


152 


15.89d=0.22 


5 




0053-016 


330 


91 


6.68±0.08 


5 




3C29 


1104 


13 


0.53±0.01 


5 


3C129 


3C129.1 





200 


20.14=b3.31 


6 




3C129 


367 


82 


6.40zb0.41 


6 


A2382 


PKS2149-158C 


300 


46 


4.08±0.10 


7 




PKS2149-158 


340 


31 


1.71±0.04 


7 


A401 


A401A 


440 


113 


7.72±0.16 


* 




A401B 


730 


74 


4.84±0.10 


* 


A2142 


A2142A 


270 


230 


53.72zb0.69 


* 


A2065 


A2065A 


1120 


48 


0.29±0.02 


* 


Ophiuchus 


OPHIB 


480 


74 


7.63±0.14 


* 



Col. 1: Cluster name; Col. 2: Source name/label; 

Col. 3: Projected distance from the cluster X-ray center; 

Col. 4: RMS of the RM distribution; 

Col. 5: Cluster X-ray surface brightness of the intracluster gas in the source location; 

Col. 6: (Trm reference: 

1 Govoni et al. (2001), 2 Bonafede et al. (2010), 3 Govoni et al. (2006), 4 Eilek & Owen (2002), 

5 Feretti et al. (1999a), 6 Taylor et al. (2001), 7 Guidetti et al. (2008), * This work. 
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